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We have observed a narrow state near 2.32 GeV/c^ in the inclusive D^tt" invariant mass distribu- 
tion from e^e~ annihilation data at energies near 10.6 GeV. The observed width is consistent with 
the experimental resolution. The small intrinsic width and the quantum numbers of the final state 
indicate that the decay violates isospin conservation. The state has natural spin-parity and the low 
mass suggests a = 0"'" assignment. The data sample corresponds to an integrated luminosity of 
91 fb~^ recorded by the BABAR detector at the PEP-II asymmetric-energy e"'"e~ storage ring. 

PACS numbers: 14.40.Lb, 13.25.Ft, 12.40.Yx 



We have found a narrow state decaying to Dfn^ at a a 91 fb ^ data sample recorded both on and off the 
mass near 2.32 GeV/c^. This result is obtained from T(4S') resonance by the BABAR detector at the PEP-II 
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asymmetric-energy e^e~ storage ring. 

Experimental information on the spectrum of the cs 
meson states is Umited. The ^S'o ground state, the 
meson, is well-estabhshed, as is the '^S'l ground state, the 
D*(2112)+. Only two other cs states have been observed 
thus far 0. The Z?si(2536)+ has been detected in its 
D*K decay mode and analysis of the D* decay angular 
distribution prefers = 0. The L»*j(2573)+ was 
discovered in its D^K^ decay mode and so has natural 
spin-parity. The assignment = 2"'" is consistent with 
the data, but is not established [3- 

The spectroscopy of cs states is simple in the limit of 
large charm-quark mass 0,0. In that limit, the total 
angular momentum j = l-\-soi the light quark, obtained 
by summing its orbital and spin angular momenta, is 
conserved. The P-wave states, all of which have posi- 
tive parity, then have j = 3/2 or j = 1/2. Combined 
with the spin of the heavy quark, the former gives total 
angular momentum J = 2 and J = 1, while the latter 
gives J = 1 and J = 0. The = 2+ and = 1+ 
members of the j = 3/2 doublet are expected to have 
smaU width 0, and are identified with the L'*j(2573)+ 
and -Dsi(2536)"'", respectively, although the latter may 
include a small admixture of the j = 1/2, = 1+ state. 
Theoretical models typically predict masses between 2.4 
and 2.6 GeV/ (? for the remaining two states j both 

of which should decay by kaon emission. They would be 
expected to have large widths 0, and hence should be 
difficult to detect. 

The experimental and theoretical status of the P-wave 
cs states thus can be summarized by stating that exper- 
iment has provided good candidates for the two states 
that theory predicts should be readily observable, but 
has no candidates for the two states that should be difh- 
cult to observe because of their large predicted widths. 

The BABAR detector is a general purpose, solenoidal, 
magnetic spectrometer, which is described in detail else- 
where 0- The detector components employed in this 
analysis are discussed briefly here. Charged particles 
are detected and their momenta measured by a combi- 
nation of a cylindrical drift chamber (DCH) and a sili- 
con vertex tracker (SVT), both operating within a 1.5-T 
solenoidal magnetic field. A ring-imaging Cherenkov de- 
tector (DIRC) is used for charged-particle identification. 
Electrons are identified and photons measured with a Csl 
electromagnetic calorimeter (EMC). 

The objective of this analysis is to investigate the 
inclusively-produced Dfn^ mass spectrum by combin- 
ing charged particles corresponding to the decay —> 
K^K~TT^ jl^ with 7r° candidates reconstructed from a 
pair of photons. Events of interest are required to contain 
at least three reconstructed tracks yielding a net charge 
of ±1 and at least two photons each of which must have 
energy greater than 100 MeV, and to have a ratio of the 
second to the zeroth Fox- Wolfram moment less than 
0.9. Charged-kaon candidates are selected based on the 



Cherenkov-photon information from the DIRC together 
with the measured energy loss in the SVT and DCH. 

A K^K^ candidate pair is combined with a third 
track that fails the kaon criteria (and so is treated as 
a pion) in a geometrical fit to a common vertex. An ac- 
ceptable K^K^TT^ candidate must have a fit probability 
greater than 0.1% and a trajectory consistent with origi- 
nating from the e"'"e~ luminous region. Background from 
K^K^ , which is evident from the corresponding 
K'^ K~ mass distribution, is removed by requiring that 
the K'^K~ mass be less than 1.84 GeV/c^. 

A candidate tt^ is formed by constraining a photon 
pair to emanate from the intersection of the K'^K^tt^ 
candidate trajectory and the beam envelope, performing 
a one-constraint fit to the 7r° mass, and requiring a fit 
probability greater than 1%. A given event may yield 
several acceptable 7r° candidates. We retain only those 
candidates for which neither photon belongs to another 
acceptable 7r° candidate. 

Finally, to reduce combinatorial background from the 
continuum and eliminate background from B-meson de- 
cay, each K^K^t:'^'!:'^ candidate must have a momen- 
tum p* in the e+e^ center-of-mass frame greater than 
2.5 GeV/c. 

The upper histogram in Fig.^a) shows the K^K^tt^ 
mass distribution for all candidates. Clear peaks cor- 
responding to D+ and mesons are seen. To re- 
duce the background further, only those candidates with 
K+K- mass within 10 MeV/c^ of the (/)(1020) mass or 
with K-'K+ mass within 50 MeV/c^ of the i?*(892) mass 
are retained; these densely populated regions in the 
Dalitz plot do not overlap. The decay products of the 
vector particles 0(1020) and i4r*(892) exhibit the ex- 
pected cos^ Oh behavior required by conservation of an- 
gular momentum, where Oh is the helicity angle. The 
signal-to-background ratio is further improved by requir- 
ing I cos 6';i I > 0.5. The lower histogram of Fig. da) 
shows the net effect of these additional selection criteria. 
The D+ signal (1.955 < m{K+K-Ti+) < 1.979 GeV/c^) 
and sideband (1.912 < m(K+K-Ti+) < 1.934 GeV/c^ 
and 1.998 < m{K+K-Ti+) < 2.020 GeV/c^) regions 
are shaded. The signal peak, consisting of ap- 
proximately 80,000 events, is centered at a mass of 
(1967.20 ±0.03) MeV/c^ (statistical error only). 

Figure ^b) shows the mass distribution for all two- 
photon combinations associated with the selected events. 
The 7r° signal (122 < m{jj) < 148 MeV/c^) and side- 
band (90 < m(77) < 110 MeV/c^ and 160 < 7^(77) < 
180 MeV/c^) regions are shaded. Candidates in the 
signal region of Fig. da) are combined with the mass- 
constrained 7r° candidates to yield the mass distribu- 
tion of Fig. n^c). A clear, narrow signal at a mass near 
2.32 GeV/c^ is seen. The shaded histogram represents 
the events in the Z?+ K^K^tt^ mass sidebands com- 
bined with the tt" candidates. In Fig. ^d) the mass dis- 
tributions result from the combination of the Df candi- 
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FIG. 1: (a) The distribution of K+K' 'n'^ mass for all can- 
didate events. Additional selection criteria, described in the 
text, have been used to produce the lower histogram, (b) The 
two- photon mass distribution from D^tt" candidate events. 

and tt" signal and sideband regions are shaded, (c) 
The Df-K^ mass distribution for candidates in the D'^ sig- 
nal (top histogram) and K'^K~n^ sideband regions (shaded 
histogram) of (a), (d) The D^77 mass distribution for signal 
Df candidates and a photon pair from the tt" signal region of 
(b) (top histogram) and the sideband regions of (b) (shaded 
histogram) . 



dates with the photon pairs from the tt*^ signal and side- 
band regions of Fig. ^Jb) (the sideband distribution is 
again shaded). In this case, all photon pairs in the signal 
region of Fig. ^b) are used. In Figs.^c) and^Jd) the 
2.32 GeV/c^ signal is absent from the sideband distribu- 
tions indicating quite clearly that the peak is associated 
with the Dfn'^ system. No other signal in the region up 
to 2.7 GeV/c^ is evident in these plots, except for a small 
i:>;(2112)+ D+TfO signal in Fig.HJc). 

In order to improve mass resolution, the nominal Df 
mass has been used to calculate the energy for the 
distributions of Fig. ^d), for the Df signal distribution 
of Fig. nfc), and for all subsequent mass distributions 
involving Df candidates. 

The Dfn^ mass distribution for p*(£)+7r") > 
3.5 GeV/c is shown in Fig. |2Ia). Similar distributions 
produced for p* values ranging from 2.5 to 4.5 GeV/c 
show the same prominent peak at the same mass value. 
The fit function drawn on Fig. |2Ia) comprises a Gaussian 
function describing the 2.32 GeV/c^ signal and a third- 
order polynomial background distribution function. The 
fit yields 1267 ± 53 candidates in the signal Gaussian 
with mass (2316.8 ±0.4) MeV/c^ and standard deviation 
(8.6 ±0.4) MeV/c^ (statistical errors only). The system- 
atic uncertainty in the mass is conservatively estimated 
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FIG. 2: The Df-jr^' mass distribution for (a) the decay 
K+K-TX+ and (b) the decay D+ ^ K+K-tt+tt^. The fits to 
the mass distributions as described in the text are indicated 
by the curves. 



to be less than 3 MeV/c^. The broad peak in Fig. [2Ia) 
centered at 2.16 GeV/c^ is due to random i:)*(2112)+7 
combinations where D*(2112)+ — > 

The signal, which we label £'*j(2317)+, is observed in 
both the (?!)7r+ and K*^K'^ decay modes of the Z?+. In 
addition, a sample of K^K^tt^tt^ decays is se- 

lected by adding tt" candidates (refit to the K'^K^n'^ 
vertex) to each K^K^n^ candidate. The purity of this 
sample is enhanced by requiring a tt*^ fit probabil- 
ity of at least 10% and selecting the K*^ , K*^ , (j), or p+ 
mass regions for the relevant two-body subsystems. Each 
resulting Df candidate is combined with a second tt'^ can- 
didate with lab momentum greater than 300 MeV/c. A 
clear £'*j(2317)+ signal is observed as shown in Fig.j^b). 
A Gaussian fit yields 273 ± 33 events with a mean of 
(2317.6±1.3) MeV/c2 and width (8.8±1.1) MeV/c^ (sta- 
tistical errors only). The mean and width are consistent 
with the values obtained for the — K'^K'^tt'^ decay 
mode. The mass distribution of the K'^ K^tt^tt^ 

sample (not shown) peaks at (1967.4 ±0.2) MeV/c^ (sta- 
tistical error only). 

We use a Monte Carlo simulation to investigate the 
possibility that the D*j(2317)"'" signal could be due to 
reflection from other charmed states. This simulation 
includes e^e~ — > cc events and all known charm states 
and decays. The generated events were processed by a 
detailed detector simulation and subjected to the same 
reconstruction and event-selection procedure as that used 
for the data. No peak is found in the 2.32 GeV/c^ 
Dfir'^ signal region. In addition, no signal peak is pro- 
duced when the and tt* identities are deliberately 
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Da (2112)"'' signal region as described in the text. The ver- 
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exchanged. 

Mass resolution estimates for the K^K^tt^tt'^ sys- 
tem are obtained directly from the data using a fit to 
the mass distribution — » if+if^vr+Tr". The mea- 
sured width from this mode is consistent with that of 



The shaded histogram corresponds to the subset of 
combinations for which either Dfj combination hes in 
the D*(2112)+ region, defined as 2.096 < m{D+"f) < 
2.128 GeV/c2. Again, no i:»*j(2317)+ signal is evident, 
thus demonstrating the absence of a 13* (2112)"'"7 decay 
mode at the present level of statistics. 

The Dfn^j mass distribution, excluding any photon 
that belongs to any tt'^ candidate, is shown as the up- 
per histogram of Fig.lS^c). The shaded histogram corre- 
sponds to the subset of combinations in which the Df^f 
mass falls in the D*(2112)"'" region. No signal is observed 
near 2.32 GeV/c^ in either case. A small peak, however, 
is visible near a mass of 2.46 GeV/c^. This mass corre- 
sponds to the overlap region of the D*(2112)"'" 
and D*j(2317)+ — *■ Dfw^ signal bands that, because of 
the smah widths of both the D*(2112)+ and £1*^(2317)+ 
mesons, produces a narrow peak in the Dfn^j mass dis- 
tribution that survives a i?*(2112)"'" selection. 

If the peak in the Dfir'^'j mass distribution of Fig.|2Ic) 
were due to the production of a narrow state with mass 
near 2.46 GeV/c^ decaying to D*(2112)+7r°, the kinemat- 
ics are such that a peak would be produced in the Dfir'^ 
mass distribution at a mass near 2.32 GeV/c^. Such 
a Dfjr'^ mass peak, however, would have a root-mean- 
square of ~ 15 MeV/c^, which is significantly larger than 
that obtained for the £'*j(2317)+ signal. In addition, 
Monte Carlo studies indicate that if the apparent signal 
at 2.46 GeV/ were due to a state that decays entirely to 
£'*(2112)+7r°, it would produce only one-sixth the signal 
we observe at 2.32 GeV/c^. 

Although we rule out the decay of a state of mass 
2.46 GeV/c^ as the sole source of the Dfir^ mass peak 



the D*j(2317)+ signal. A simulation of the D*j (23 17)+ corresponding to the D*j(2317)+, such a state may 

be produced in addition to the -D*j(2317)+. How- 



decay to K'^ K~Tr~^TT'^ yields a similar mass resolution af- 
ter event reconstruction and selection criteria have been 
satisfied. We conclude that the intrinsic width of the 
D*j(2317)+ is small (T < 10 McV). 

The cosOh distribution of the i:';j(2317)+ decay with 
respect to its direction in the e+e~ center-of-mass frame 
has been investigated. The efhciency-corrected distribu- 
tion is consistent with being flat, as expected for a spin- 
zero particle, or for a particle of higher spin that is pro- 
duced unpolarized. 

We have also performed a search for the decay 
13:^(2317)+ ^ L»+7. Shown in Fig. Ela) is the D+7 
mass distribution obtained by combining a candi- 
date in the signal region of Fig. QJa) with a photon with 
an energy of at least 150 MeV that does not belong to a 
77 combination in the signal region of Fig. ^b). The re- 
quirement that the p* of the Dfj system be greater than 
3.5 GeV/c is also imposed. There is a clear Z)*(2112)+ 
signal, but no indication of £'*j(2317)+ production. 

The -0^77 mass distribution for p* (13+77) > 
3.5 GeV/c, excluding any photon that belongs to the 7r° 
signal region of Fig.^b), is shown as the upper histogram 
of Fig. ISIb). No signal is observed near 2.32 GeV/c^. 



ever, the complexity of the overlapping kinematics of the 
£i:(2112)+ ^ £1+7 and i:»;j(2317)+ ^ D+Tr° decays re- 
quires more detailed study, currently underway, in order 
to arrive at a definitive conclusion. 

The decay of any cs state to -D+tt" violates isospin 
conservation, thus guaranteeing a small width. It is pos- 
sible that the decay proceeds via 77 — tt" mixing, as dis- 
cussed by Cho and Wise 0. For a parity-conserving 
decay only a spin-parity assignment in the natural se- 
ries {0+, 1~, 2+, . . . } is allowed. The low mass compared 
to those of the Dsi(2536)+ and the 13*^(2573)+ favors 
= 0+. In this case, decay to Dfj is excluded. How- 
ever, decay of the 13*^(2317)+ to D;(2112)+7 is allowed 
and might compete with decay by pion emission. The 
shaded mass distribution of Fig. Efb) suggests that this 
mode is absent, at least at the present level of statistics. 
This may simply indicate that decay by pion emission is 
favored over radiative decay. 

Further studies are under way. If, however, the tenta- 
tive = 0+ assignment is confirmed, the low mass, 
small width, and decay mode of the Z3*j (23 17)+ are 
quite different from those predicted by potential mod- 
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elsilli. 

In summary, in 91 fb ^ of data collected by the 
BABAR experiment we have observed a narrow state in 
the inclusive Dfir'^ mass distribution near 2.32 GeV/c^. 
We find no evidence for the decay of this state to 
£'*(2112)+7, or Z3+77. Since a cs meson of this 
mass contradicts current models of charm meson spec- 
troscopy (fll, |3, S| J either these models need modification 
or the observed state is of a different type altogether, 
such as a four-quark state. 
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